Photonic band structure of ZnO photonic crystal slab laser 
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We recently reported on the first realization of ultraviolet photonic crystal laser based on zinc 
oxide [Appl. Phys. Lett. 85, 3657 (2004)]. Here we present the details of structural design and 
its optimization. We develop a computational super-cell technique, that allows a straightforward 
calculation of the photonic band structure of ZnO photonic crystal slab on sapphire substrate. 
We find that despite of small index contrast between the substrate and the photonic layer, the 
low order eigenmodes have predominantly transverse-electric (TE) or transverse-magnetic (TM) 
polarization. Because emission from ZnO thin film shows strong TE preference, we are able to limit 
our consideration to TE bands, spectrum of which can possess a complete photonic band gap with 
an appropriate choice of structure parameters. We demonstrate that the geometry of the system 
may be optimized so that a sizable band gap is achieved. 
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PACS numbers: 42.55.Tv,78.66.Hf,42.55.Px 



I. INTRODUCTION 

A photonic crystal slab (PhCS) is a layer of dielectric 
where refractive index is periodically modulated within 
the plane of the slabQ, 0, H- Such devices have at- 
tracted much attention because of their potential ap- 
lications to various optoelectronic devices and circuits 

' H IS II S H E3 Phcs ma y lead t0 a complete 
photonic band gap for the guided modes when the wave 
propagation in all in-plane directions is forbidden due 
to Bragg interference. Moreover, its planar geometry 
makes it easy to incorporate on a chip, as well as of- 
fers a possibility of low-threshold H E Jll H H E3 
electrically-driven compact laser source lfij. Defect cavi- 
ties in these photonic devices can have high quality fac- 
tor and small modal volume 0, 0, 0, Q|. So far the 
experimental efforts mainly concentrated on PhCS made 
of III-V semiconductors. [E [H [H [H Q 13 E3 They 
operate in the infrared (IR) communication frequencies 
(see also 0), where the structural feature size is on the 
scale of 0.25/im. 

There is technological and commercial demand for 
compact and integrable laser sources in near the ultra- 
violet (UV) range of optical spectrum. Recently, we re- 
ported the first realization of UV PhCS laser that oper- 
ates at room temperature |2 lj . This required overcoming 
several experimental and design challenges. 

First of all, because of the shorter wavelength, PhCS 
requires smaller structural features with sub lOOnm size. 
The fabrication of such miniature structures is technolog- 
ically challenging for commonly used wide band gap semi- 
conductors such as GaN and ZnO. Compared with the 
other wide band gap materials ZnO, has the advantage 
of large exciton binding energy (~ 60 meV), that allows 
efficient excitonic emission even at room temperat ure 22) . 
Using Focused Ion Beam (FIB) etching technique |2l| we 
were able to achieve more than a four-fold reduction, 
compared to IR PhCS, of the feature size in ZnO (see 
also H|). 



Secondly, high quality ZnO films are usually grown 
on the lattice matched sapphire substrate |24j. that re- 
sults in relatively low refractive index contrast between 
PhCS and substrate. Furthermore, as we detail below, 
the substrate-PhCS-air rather then air-PhCS-air geom- 
etry in vertical direction disallows the classification of 
the photonic modes into symmetric and antisymmetric 
with respect to reflection in the mid-PhCS plane 0,0113. 
Separation of the guided modes according to their sym- 
metry in IR PhCS is usually ensured due to "under- 
cutting" , i.e. removing the material directly below the 
photonic layer 1!% 1 151 llfij . The separation can also be 
approximatelyj3,|26j done for systems with a large index 
contrast between the substrate and PhCS. Neither option 
is readily available for ZnO-based structures, due to its 
material properties. At first glance, this makes photonic 
band gap (PBG) unattainable in our system. The pur- 
pose of this paper is to design ZnO-based PhCS, which 
exhibits PBG effects. 

In Section II we introduce a computational super-cell 
technique that allows the calculation of photonic band 
structure of PhCS on a dielectric substrate. Then we 
present the detailed analysis of symmetry and polariza- 
tion of the photonic modes belonging to different bands. 
We show that low-order photonic bands can still be ap- 
proximately classified as TE and TM. In Section III we 
make use of the experimental fact that PhCS is made 
of single-crystalline hexagonal ZnO thin film, which is 
grown along c- crystallographic axis [2^| . Since the polar- 
ization of exciton emission in the ZnO thin film is pre- 
dominantly perpendicular to the c-axis 27, 28], the emis- 
sion goes preferably into TE-polarized photonic modes. 
Coincidentally, the air-holes-in-a-dielectric PhCS config- 
uration can possess complete PBG (with certain param- 
eters) for namely TE polarized modes. Therefore, we 
consider only TE-polarized photonic modes in our band 
structure calculation. Finally, in Section IV we systemat- 
ically study the position and the size of PBG as a function 
of geometrical parameters: air-hole radius r, slab thick- 
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FIG. 1: Schematic diagram of the vertical super-cell used 
in photonic band structure simulation. Cylindrical air holes 
of radius r are arranged in a hexagonal pattern with lattice 
constant a. The holes extend throughout the slab of thickness 
tphcs an d may penetrate into the substrate to the depth of 
t e . The thicknesses t a and t s need to be chosen sufficiently 
large (see text for discussion). The cell is symmetric with 
respect to reflection in E sc plane. For convenience we also 
introduce a middle plane of the slab, E. 

ness tphcs and etch depth t e as illustrated in Fig. We 
conclude our paper with a discussion section. 



II. SUPER-CELL PHOTONIC BAND 
STRUCTURE CALCULATION 

Although "Bragg reflection" plus "index guiding" de- 
scription of the electromagnetic field in PhCS can eluci- 
date many physical phenomena that occur in such sys- 
tem, it is naive because both in-plane and vertical propa- 
gation or confinement are inseparable. In order to obtain 
the photonic modes of PhCS one needs to solve the prob- 
lem self-consistently within the same framework. This 
can be done with a plane- wave expansion method [2^] that 
de facto became a standard in the field. 

PhCS designed for IR applications have an inherently 
large size of the structural units owing to relatively long 
wavelengths involved. In these systems it is possible to 



fabricate PhCS cladded with air above and below (see e.g. 
[MEjEI)- This design is advantageous for a number of 
reasons. First, this is a symmetric structure - there ex- 
ists a reflection symmetry with respect to the plane that 
goes through the middle of PhCS. This symmetry allows 
two orthogonal classes of eigenfunctions of Maxwell equa- 
tions: symmetric and antisymmetric Q. Super-cell for 
the band-structure calculation can be easily constructed 
by repeating ...-air-PhCS-air-PhCS-... indefinitely in the 
vertical direction 3] . It was shown that low-lying (E-field) 
symmetric modes are predominantly TE polarized, while 
antisymmetric ones are TM polarized. Due to superior 
guiding properties of TE modes and the fact that self- 
supporting membrane percolated by air holes (unlike col- 
lection of dielectric cylinders in air) exhibits PBG for the 
same polarization, this geometry is widely used in prac- 
tice. The super-cell photonic band structure calculations 
also give unphysical modes that do not correspond to 
the guided modes of PhCS (confined to the slab). For- 
tunately, physical (guided) modes can be separated us- 
ing the concept of light-cone as follows. In the result of 
photonic band structure calculation one obtains the fre- 
quency u> for a given in-plane wavevector k| | . This EM 
mode (concentrated inside the photonic layer) is coupled 
only to the modes outside the layer with the same (w,k||). 

If u> < cfc||, then k z = ^Ju) 2 /c 2 — is imaginary outside 

PhCS. Thus the light intensity decays exponentially in 
the z-direction away from the layer. This describes a 
guided mode. In the opposite case, k z is real and light 
can escape to the infinity Such a mode is a leaky mode, 
which cannot be accurately described in the super-cell 
calculations. A continuum of modes with to > cfc|| form 
a light cone (a cone in u> — k space). Therefore, PBG 
obtained in all planar geometries are, strictly speaking, 
photonic band gaps for guided modes. 

In our UV PhCS laser [53 we used ZnO because of 
its wide electronic band gap and robustness of the op- 
tical properties with respect to damage caused by FIB. 
There are two consequences relevant for the band struc- 
ture calculation. First, in order to produce high qual- 
ity films, ZnO was grown epitaxially on sapphire sub- 
strate with large refractive index n s ~ 1.78. Second of 
all, our choice of the materials and miniature structure 
dimensions, do not currently allow for removal of the un- 
derlying layer of sapphire — "undercutting" is impossi- 
ble. These experimental limitations lead to severe re- 
strictions, (i) Presence of the sapphire substrate lifts the 
reflection symmetry This results in mixing between 
TE-like and TM-like classes of electromagnetic modes. 
One may expect that this effect should be large in view 
of relatively small index contrast between sapphire and 
ZnO {nznO = 2.35). Below we will show that this ex- 
pectation is not completely borne out. (ii) The most 
stringent restriction is imposed by the presence of the 
second light cone - that reflects the possibility of radia- 
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FIG. 2: Reflection symmetry allows for solutions symmetric 
and antisymmetric with regard to reflection in E ac . This arti- 
ficial symmetry helps to separate the guided modes from the 
leaky ones, see text for details. 

tive escape into the sapphire substrate. Indeed, although 
light leakage into the air occurs for u> > efei i , leakage into 
the substrate starts already for oj > (c/n s )fcn. Avoiding 
this problem requires the use of high filling fraction in 
ZnO PhCS (/ > 50%), that in turn drives down the size 
of the features to be produced in the experiment. 

A band structure calculation technique for the asym- 
metric geometry like ours was proposed in Ref. Q . The 
authors of Ref. 0] proposed to use a super-cell where 
substrate-PhCS-air-substrate-PhCS-air-.. is periodically 
repeated in z-direction. Special care needs to be exercised 
in order to avoid unphysical solutions that correspond 
to the guided modes in the substrate layers bounded 
by air and a PhCS. In our simulation we constructed 
a different super-cell by combining the structure and its 
mirror image (shown in Fig. ^) for the photonic band 
calculation 29] . Our choice of the unit cell is motivated 
by the following. Due to particular arrangement of the 
layers in our super-cell, there exists an artificial symme- 
try - reflection plane S sc in Fig. ^ Therefore, there 
is an artificial classification of the modes into symmetric 
and antisymmetric, as it is schematically shown in Fig. 
12 This property provides a way to control the precision 
of the band structure calculation. Indeed, only for the 
PhCS guided modes (with eigenfrequencies outside of air 
and substrate light cones), and only when thicknesses of 
air (t a ) and substrate (t s ) layers are chosen to be suffi- 
ciently large, the eigen-frequency of S sc -symmetric and 
E sc -antisymmetric modes are nearly degenerate. This 
can be seen from Fig. there should be no difference 
in eigenenergy between the functions when the value of 
the eigenfunction is close to zero at (i) E sc plane and (ii) 
the top/bottom of the super-cell. The second condition 
arises because of the periodicity in the vertical direction 
— the top/bottom of the super-cell is also a middle plane 



that separates two closest PhCS layers. 

Furthermore, one may think that the doubling of the 
1600 1 i 
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FIG. 3: Experimentally measured emission intensity from 
600nm ZnO film as a function of polarization angle. The 
maximum intensity at about 25° corresponds to E-field po- 
larization parallel to the plane of the film (TE polarization). 
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FIG. 4: Calculated band structure of ZnO photonic crystal 
slab air cylinder radius r/a — 0.25, and slab thickness t/a — 
1.4. The refractive indexes for ZnO and sapphire are 2.35 and 
1.78. Thin and thick lines represent TM and TE polarized 
modes respectively. 



cell size should increase the computation time. However, 
at least for a PhCS with a hexagonal symmetry, S sc con- 
tains a center of inversion. The latter allows to perform 
the photonic band structure calculation assuming inver- 
sion symmetry that shortens the calculation time (of one 
symmetry) by a half. In all simulations reported below, 
we checked the consistency of the results with the pro- 
cedure outlined above, then studied the band structure 
of the modes outside of the substrate (more restrictive) 
light cone. 
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FIG. 5: Integrated according to Eq.Q x, y and z components 
of the electric, E, and magnetic, H, field. Solid/dashed lines 
represent real/imaginary part of the field. First /second two 
columns correspond to the first/second mode of Fig. 2] taken 
for fc|| a,t K symmetry point. The same normalization is used 
for all field components, so that their values can be compared. 
Mode M1/M2 shows a pronounced TE/TM polarization. 



M. 



-1000 
-2000 

-3000 
-4000 
-5000 

-60og 

-500 

-1000 

-1500 

-2000 

-2500 

-3000 
0.04 
0.02 


-0.02 
-0.04 
-0.06 
-0.08 
-0.1 
-0.12 
-0.14 
Air Sbtr 



X \ 




< H y> 




<H z >^- 

r \ 

\ 

\ 

1 

1 

f 
1 

1 ' 


1 

1 

1 

1 



PhCS 



Air 



III. PHOTONIC MODE CLASSIFICATION: 
POLARIZATION 



We measured the polarization of ZnO emission inten- 
sity from a 600nm-thick ZnO film on sapphire substrate, 
c-axis of ZnO is perpendicular to the surface of the film. 
ZnO is optically excited by the third harmonic of a mode- 
locked Nd:YAG laser (355nm, 10Hz, 20ps). A cylindrical 
lens (/ = 100mm) is used to focus the pump beam to a 
strip normal to the edge of the sample. The emission was 
collected at 1/iJ pump pulse energy and the pump area 
was about 2mm x 50fj,m. Fig. clearly demonstrates 
that the emission is strongly polarized with the elec- 
tric field parallel to the film (TE polarization) p H E ^ . 
Because the emission in the ZnO thin film is predomi- 



nantly TE polarized we need to consider consider only 
the TE-polarized photonic modes. This motivates a de- 
tailed analysis of the polarization of the modes in the 
experimentally relevant PhCS. 

In Fig. 0]we show an example of calculated photonic 
band structure. Air light-cone boundary, uj — ckn, and 
substrate light cone boundary, uj — (c/n s )k\\, are shown 
with dash-dot and dashed lines respectively. The pres- 
ence of sapphire cone limits any bandgap to a/A < 0.325. 
For PhCS laser operation we need to overlap a photonic 
band gap with gain spectrum of ZnO, 380nm < A < 
AOOnrn. Consequently, the lattice constant in the hexag- 
onal air-hole pattern of PhCS has an upper bound of 
a < 130nm. To fabricate such small features we used the 
FIB etching^. 
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Since the presence of the substrate removes the sym- 
metry with respect to reflection in PhCS middle plane 
(E in Fig. QJ, one cannot separate the modes into two 
independent classes and consider them separately. This 
separability was crucial for obtaining sizable PBG in a 
photonic membrane (or strictly two-dimensional struc- 
tures) . For comparison, we calculated the band structure 
of the system with the same parameters as in Fig. SI but 
sapphire substrate was replaced by air. The obtained dis- 
persions looked qualitatively similar to the ones in Fig. 

so that it was possible to make band-to-band corre- 
spondence. This has been used previously to justify the 
existence of PBG in the spectrum of TE modes Q ■ In our 
case the refractive index of the substrate is quite large, 
moreover it may even become comparable to the effective 
index of PhCS as the patterning of ZnO film with orig- 
inal refractive index nznO = 2.35 will reduce it to just 
a little above that of the sapphire. These considerations 
warrant an in-depth analysis of the mode polarizations. 

In order to determine the degree of polarization of the 
electromagnetic eigenmodes of PhCS, we studied 



(F(z)) 



F(x,y,z) dx dy, 



(1) 



where F represents x, y or z component of electric or 
magnetic field. Fig. [3] shows these quantities for the 
first two bands taken at M-symmetry point of Fig. 01 
It is clearly seen that the first (the lowest frequency) 
mode Mi is strongly TE polarized: E x , E y and H z 
are more than one order of magnitude larger the other 
components. Furthermore, as discussed above dominant 
components are nearly symmetric with respect to middle 
plane of PhCS, whereas H x , H y and E z are close to an- 
tisymmetric. The deviation from the perfect symmetry 
is entirely due to the presence of the sapphire substrate. 
Comparing the second mode M% to the previous case one 
notices that the polarization of the field has changed to 
the opposite — M 2 is predominantly TM polarized. 

An increase in the order of the band does lead to en- 
hancement of polarization mixing. However as we will 
see in more detail below, the substrate light-cone con- 
dition restricts the involvement of higher-order bands. 
We checked the polarization of the lowest six bands in 
several k\\ points and found that in all relevant cases 
certain symmetry, TE or TM, can be assigned to each 
band. Therefore, despite of highly asymmetric claddings 
(sapphire versus air) and relatively small refractive index 
contrast between PhCS (n e ff ~ 2.1) and the substrate 
(n s = 1.78), the mixing of the polarizations in the low- 
order photonic bands is limited. From now on we will 
ignore TM polarized bands. In the next section we will 
systematically study the effect of the structural param- 
eters on the width of PBG in the spectrum of TE (see 
thick solid lines in Fig. 0} modes. 
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FIG. 6: The position and relative size of photonic bandgap 
for PhCS with variable structural parameters, (a) tphcs/ a = 
1.4, t e = 0; (b) r/a = 0.25, t e = 0; (c) r/a = 0.25, 
tphes/o- = 1.4. riznO = 2.35 and n s = 1.78 were used in 
the calculation. The left vertical axes should be used with 
solid lines, which represent the edges of PBG. Dashed lines 
and the right vertical axes show the relative size of the gap. 



IV PHOTONIC BAND GAP OPTIMIZATION 

We intend to optimize the following three parameters 
that can be controlled in the FIB etching 21]: (i) PhCS 
thickness, tphes, (h) air hole radius, r (related to the 
filling fraction /) and etch depth, t e . Due to linear nature 
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of Maxwell equations, we normalize all parameters to the 
value of lattice constant, a. Our goal is to maximize the 
relative (normalized to the center frequency) PBG width. 
Once this is achieved, we can overlap PhCS gap with 
the emission spectrum of ZnO by choosing appropriate 
a. From the experimental prospective it is difficult to 
reproduce extremely small structural features, therefore 
we prefer to have PBG at high values of a/ A if possible. 

Fig. shows the results of our simulation. In panel 
(a) of the figure, we plot the position and relative width 
of PhCS as we vary r/a within the range 0.18 — 0.30, 
and keep tphcs/ a — 1-4 and t e = 0.0. This corresponds 
to a change in the pattern filling fraction from / = 0.88 
to / = 0.67. One can clearly see that there exists a 
PBG maximum at r/a ~ 0.24 (or / ~ 0.79). This ef- 
fect has a clear physical interpretation: at large filling 
fractions (small r/a) the index contrast within PhCS is 
small, therefore band splitting, as well as PBG, is small. 
At small filling fractions, the wave guiding becomes poor, 
furthermore the decrease of the effective refractive index 
of the photonic layer leads to a shift of PBG as a whole 
to higher frequency. Eventually, at r/a ~ 0.26 the top of 
the gap is set not by the higher order guided mode but 
by a/A ~ 0.324 at the light-cone when k\\ is at M-point 
(see Fig. 0J. In reality, above this frequency the guided 
modes will always be coupled to some leaky modes inside 
the light-cone. This phenomenon results in a plateau of 
the upper bandgap edge for r/a = 0.26 — 0.30. Conse- 
quently, PBG disappears at r/a > 0.30. For the air holes 
with such a large radius, the effective refractive index of 
PhCS n e ff ~ 1.87 approaches that of the substrate. 

Fig. EJb) presents the effect of the photonic layer thick- 
ness, tphcsi on PBG. With the increase oit Ph cs we also 
observe a maximum. For samples thinner than the corre- 
spondent wavelength, the eigenmodes of the system sub- 
stantially "spill" into the air and (mostly) the substrate. 
Consequently, the effective index of PhCS becomes com- 
parable to n s so that PBG lies at higher frequencies. The 
light-cone leads to the plateau of the upper bandgap edge 
as in the case of large r/a in Fig. HJa). For very thick 
photonic layers, tphcs/ a > 1-5, the gap closes again. 
This occurs due to a stronger dependence of upper band 
edge frequency on the effective refractive index of PhCS. 
Indeed, as n e f f is increased the band forming the upper 
edge of PBG develops the kinks close to the light-cone, 
that leads to the reduction of PBG as seen in Fig. [Bfb). 
As a result of these dependences, the size of the gap can 
be maximized at tphcs/a ~ 1.45. 

The possibility of improving guiding properties of the 
photonic layer by extending the air holes into the sub- 
strate was discussed in Ref. [3(|- Indeed, removing a 
part of the substrate material lowers its refractive index 
and achieves a partial effect of "undercutting" . However, 
the fact that PBG in our case lies close to the light-cone 
has a significant effect. An increase in t e results in the 
decrease of the effective refractive index experienced by 




r/a 




FIG. 7: The degree of polarization of the modes calculated 
for fen in M (solid lines) and K (dashed lines) k-points. The 
modes belonging to the lowest four bands are considered. 
Cases (a)-(c) correspond to those of Fig. HJ (...) denotes 
an integration over x and y coordinates. Discontinuities in 
the curves occur when the order of bands changes. 

a guided mode, because it extends into the substrate. 
Therefore the eigenfrequency of the mode increases and 
may enter the light-cone 3 lj . This is exactly what hap- 
pens in our system as shown in Fig. Et- As t e is increased 
the upper edge of PBG becomes defined by the position 
of the substrate light-cone at M k-point (Fig. 0J, whereas 
the mode that defines the lower edge of PBG is more con- 
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fined to PhCS and, therefore, is less sensitive to the effect 
of t e . As a result, Fig. EIc) shows a weak dependence on 
the etch depth. 

As we discussed in the previous section our calculation 
of the photonic band gap relies on the possibility of sep- 
arating the photonic bands based on their polarization. 
To check the consistency of this assumption we calculated 
the ratios between in-plane and normal components of 
the electric fields for various structural parameters. In 
Fig. [7| we show the polarization ratios calculated for the 
lowest four bands at two ku vectors (M and K k-points). 
As the structural parameters are varied, the polarization 
of the modes is preserved. The most noticeable change 
occurs with the increase of r/a in Fig. Eta). This is due 
to significant modification of the band structure that oc- 
curs in this case. Notwithstanding, the modes still re- 
main strongly polarized, justifying our methodology. 



DISCUSSION AND CONCLUSION 



designed pattern with 

a ~ 123ron; tphcs — 180rtm; t e ~ Onm; r ~ 30nm. (3) 

This has been achieved with the FIB etching technique. 
In order to avoid an additional damage to PhCS we did 
not attempt to remove the substrate, in other words, ex- 
perimentally we preferred t e = 0. The maximum rela- 
tive width of PBG that can be achieved via optimization 
is 5%. This is significantly smaller than what is b eing 
used in the membrane (air-bridge) geometry^, li p . 1 1 fi ll 
or in the case of low refractive index substrate|7j, l2fij . 
Narrow gap makes it difficult to align an intentionally 
introduced defect mode inside PBG. However, there are 
always some defect modes with frequencies inside PBG 
formed by small disorder unintentionally introduced in 
fabrication process. These are the lasing modes observed 
in our experiment. 

This work was supported by the National Science 
Foundation under the grant no. ECS-0244457. 



In order to simulate the realistic structures to be 
realized in the experiment |2lj. we modified super-cell 
technique^ within plane- wave expansion method for the 
photonic band structure calculation^. We constructed 
a symmetric super-cell with the doubled size, Fig. ^ The 
induced degeneracy of eigenvalues provides a convenient 
way to verify the self-consistency of the calculation. 

Substrate with high refractive index was expected to 
significantly mix the polarization of eigenmodes of PhCS. 
However, our calculations demonstrate that although the 
substrate induces asymmetry of wavefunctions (Fig. [SJ, 
they remain strongly TM- or TE-polarized for low-order 
bands. High filling fraction air-hole-in-ZnO-matrix ge- 
ometry can possess a complete PBG for TE bands, mean- 
while, ZnO film with c-axis along the growth direction 
emits mainly into TE-polarized modes. This enabled us 
to build a UV PhCS laser [i^- In the photonic band 
structure calculation we find the optimum set of param- 
eters for maximum PBG: 

1-0.24 t£0£l~iA5 -~0.0 (2) 
a a a 

These parameters are significantly different from the typ- 
ical parameters for IR PhCS. This is due to several fac- 
tors: (i) presence of the sapphire substrate brakes vertical 
symmetry of PhCS; (ii) the refractive index contrast of 
ZnO / sapphire is lower than that of InP / air usually used 
in IR PhCS; (iii) to preserve guiding in the photonic layer 
the filling fraction and the thickness of ZnO PhCS needed 
to be significantly increased. This in turns explains the 
relatively small PBG of about 5% that can be obtained 
in the ZnO PhCS on sapphire. 

Overlapping the photonic bandgap with the emission 
spectrum of zinc oxide required precise control of the 
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